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^ ! We address the question of whether anisotropic superconductivity is compatible 

s, 

, with the evidently weak sensitivity of the critical temperature Tc to sample quality in 

C , 

Q , the high- Tc copper oxides. We examine this issue quantitatively by solving the strong- 

O 

coupling Eliashberg equations numerically as well as analytically for s-wave impurity 
scattering within the second Born approximation. For pairing interactions with a 
characteristically low energy scale, we find an approximately universal dependence of 
the d-wave superconducting transition temperature on the planar residual resistivity 
which is independent of the details of the microscopic pairing. These results, in 
conjunction with future systematic experiments, should help elucidate the symmetry 
of the order parameter in the cuprates. 
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A growing body of experimental evidence has been interpreted as supporting an 
anisotropic pairing state in the high-temperature superconductors There are indica- 
tions, however, that the measured values of the superconducting critical temperature Tc do 
not depend strongly on sample quality. Because both magnetic and non-magnetic impuri- 
ties are pair-breaking in anisotropic superconductors 0, this latter statement appears to be 
incompatible with anisotropic pairing. 

Two arguments can be offered to explain this apparent contradiction, but neither is 
quantitative enough to settle the issue. First, the coherence lengths ^ 10 A in the 
cuprates, while the mean free paths I ~ 100 - 200 A [^]. Since one expects to be 
affected by impurities only when ^ ~ /, the small ratio ~ 0.1 suggests that Tc should 
not be sensitive to defects. However, for magnetic impurities in s-wave superconductors, 
superconductivity is completely destroyed for ^/l = 0.12 - 0.17 0], implying that in the 
cuprates may not be small enough to avoid an impurity- induced reduction of Tc. Second, one 
observes significant inelastic scattering in these materials which could act to mask impurity 
effects. Inelastic scattering would only be able to screen the defects effectively if the inelastic 
mean free path were much shorter than the impurity mean free path. Estimates of these 
quantities indicate that they are comparable in high-quality samples. 

In this paper, we examine quantitatively the effect of impurities in the cuprates by gener- 
alizing the Abrikosov-Gor'kov scaling law |^ to anisotropic strong-coupling superconductors. 
We present model-independent predictions of Tc as a function of planar residual resistivity 
due to non-magnetic impurities for a 90 K d-wave superconductor and compare these results 
to the response of an s-wave superconductor. We also check the validity of this analytical 
result by calculating numerically the suppression of Tc induced by structureless impurities 
treated within the second Born approximation and in strong-coupling Eliashberg theory. 
Finally, we discuss how these predictions may be tested experimentally. Other authors have 
examined this question using analytical as well as numerical techniques . The goal 

of the present work, however, is to quantify the expected Tc suppression for d-wave models of 
the high temperature superconductors, so that the d-wave hypothesis can be more directly 
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tested. 

We compute T^. in the presence of impurities from the standard mean field formahsm 
in which the electron self-energy is solved self-consistently from the single-exchange graph. 
This approach is justified if Migdal's theorem |T0[ applies, which is the case in conven- 
tional, phonon-mediated superconductors |TT]], but it is not clear whether a similar results 
holds in the high-temperature superconductors Nonetheless, we follow other authors 



i6|j7|J^,p!^, |T^ , |T5[] and assume that this result holds in what follows. 

In most of our calculations, we employ the standard set of approximations to this mean- 
field theory in order to simplify the resulting Eliashberg equations These simplified 



equations are discussed in more detail in Ref. [0. The central approximation in this ap- 
proach is to limit the wavevectors of the electron self-energy and pairing potential to the 
Fermi surface, so we refer to the solution of these approximate equations as the Fermi- 
surface-restricted solution. For comparison, we have also solved the Eliashberg equations 
without these approximations on two-dimensional lattices of small size (64 x 64, 32 x 32) 
using recently developed fast Fourier transform techniques |^,0,|T^ . We refer to the solution 
of these equations as the exact solution. 

Of the models of d-wave superconductors available , we employ a spin- 

fiuctuation-mediated pairing interaction with a spectrum constrained by neutron scattering 
in YBa2Cu306.7 [0. Solving the Fermi-surface-restricted equations with the pairing poten- 
tial and band structure discussed in Ref. [0 gives us the critical temperature as a function 



of the bare impurity scattering rate t^^. (Throughout this paper, we set h = ks = 1). We 
compute Tc vs. T~^p for a variety of critical temperatures in the absence of impurities T^o and 
in both the weak- (inelastic scattering ignored) and strong- (inelastic scattering included) 
coupling cases. Although we vary the T^o's in this paper, constraints on the electron-spin 
fluctuation coupling constant imposed by ac conductivity require that T^o = 7.2 K in this 
model |T^. For comparison, we also compute the impurity- induced Tc suppression from 
the exact Eliashberg equations for the Monthoux-Pines spin-fluctuation model PJTB[] with a 
coupling constant chosen so that Tco = 100 K. 



We will show that these numerical results can be approximately collapsed to a universal 
curve. One can deduce that the form of this curve for d-wave superconductors is 

Tr,\ . / 1 aTco \ . / 1 



where ip{z) is the digamma function, a = 1/{2{1 + Xz)TmpTco) is the strong-coupling 
pair-breaking parameter, and Az = 1 - {Im[I]''°'^°"-'(k,iLJo)])/LJo. In this last equation, 
Y^boson^^^ ztuo) refers to the self-energy of the pairing boson at the lowest fermionic Matsubara 
frequency uq = nTco, and the angle brackets denote an average over the Fermi surface. We 
note that, for an s-wave superconductor with magnetic impurities, Eq. (|ip would still apply, 
but without the factor of two in a. We note that individual aspects of this formulation have 
appeared in the works of other authors [0,^,0; in particular, this form of the pair-breaking 
parameter can be found in Ref. [§. 

In order to make contact with experiment, we represent the pair-breaking parameter a 
in terms of the planar residual resistivity po- From the Kubo formula under the standard 
assumptions that the pairing boson energy is small compared to the electronic energies and 



that vertex corrections are not important [|18[, the real part of the low frequency residual 
electrical conductivity is given by 

ReaM = ^— (2) 

In this equation, u*f = oj^J (1 + A) is the renormalized plasma frequency which is measured 
experimentally, 

1 1 1 



'^imp 1 + A Timp 



(3) 



is the renormalized scattering rate due to impurities, and A = — (9ReS "(k, 

uj=0 

is the mass renormalization parameter due to the pairing bosons. At zero temperature, A 
is equal to the inelastic scattering parameter from the Eliashberg equations Xz- Since the 
characteristic energy of the pairing boson is low compared to other electronic energy scales 
(i.e., of the order of T^), we find that A = A^ to within 5-10 % at T^. Combining this result 
with Eq. @) and the definition of the pair-breaking parameter a, we see that 



An 

Po = -^2Tco«. (4) 

We emphasize that Eq. allows us to predict the suppression as a function of planar 
residual resistivity from the experimentally observed plasma frequency uj*f and the critical 
temperature Tco independently of any microscopic model. 

Having set up the formalism for our calculations, we will now discuss the results. Fig. 
1 shows Tc/Tco vs. scattering rate from non-magnetic impurities for a selection of Tco's in 



the model d-wave superconductor of Ref. [0. Since Tco is inversely proportional to the 
coherence length ^, the different curves in each figure correspond to different ^'s. Both 
weak- (Fig. 1(b)) and strong- (Fig. 1(a)) coupling calculations show that, as T^q increases 
(^ decreases), the superconductor becomes less sensitive to the presence of non- magnetic 
impurities. In addition, by comparing curves in Figs. 1(a) and 1(b) with the same Tcq, 
we see that the strong-coupling curves (which include the effects of inelastic scattering) are 
less sensitive to impurities than the weak-coupling results |]T9[. These trends conform to the 
qualitative expectations discussed in the introduction. 

To confirm the validity of our analytical results for the scaling function Eq. (|l]), we 
plot in Fig. 2 the numerical data in Fig. 1 in terms of the scaled variables Tc/Tco and a. 
As can be seen from the figure, all of the scaled curves cluster near the scaling function. 
Note that the data in Fig. 2 include both weak- and strong-coupling results; moreover, even 
the curve from Fig. 1(a) with nearly the maximum achievable T^o = 50 K falls near the 
scaling curve. In the inset to Fig. 2, we display the impurity-induced suppression of a 100 K 
d-wave superconductor in the Monthoux-Pines model computed from the exact Eliashberg 
equations. When plotted in terms of T^/Tco and a, this model also produces a Tc suppression 
which is close to the scaling function Eq. (|1|). 

Having established the validity of the scaling law, we use this relation to predict the 
response of a 90 K superconductor to non-magnetic impurities. In Fig. 3, the shaded region 
corresponds to the d-wave Tc computed from Eqs. (0) and (H) for plasma frequencies u*i 
ranging from 1.1 to 1.4 eV as a function of the planar residual resistivity po- These plasma 



5 



frequencies are chosen to reflect tlie range of experimental uncertainty in uj*i in YBa2Cu307 
pO| . For comparison, Fig. 3 also shows the expected response of an s-wave superconductor 
to non- magnetic impurites (dashed line); in accordance with Anderson's theorem ||21]], these 
impurities have no effect on Tc. 

In addition to the uncertainty in the plasma frequency, the accuracy of the prediction in 
Fig. 3 is affected by the applicability of the scaling law and the assumption of structureless 
impurities. From Fig. 2, one can see a systematic trend away from the scaling curve as 
Tco increases. For a 90 K superconductor, this error amounts to roughly a 20 % correction 
to the horizontal scale in Fig. 3, which is about the same magnitude as the uncertainty 
in the plasma frequency. The effect of ignoring the wavevector dependence of the impurity 
scattering matrix element is more difficult to estimate. In extreme cases, impurities with 
wavevector structure could alter the scale on the horizontal axis by a factor of two . Even 
considering these caveats, if the superconducting order parameter has d-wave symmetry, 
then the prediction in Fig. 3 should give the correct scale for the residual resistivity at 
which significant depression of the critical temperature will occur in YBa2Cu307 . 

Experimentally, one can estimate the residual resistivity po by extrapolating the mea- 
sured linear resistivity vs. temperature for temperatures T greater than to T = 0. It is 
found that this extrapolated planar po is about 20 pf2-cm. in high-quality twinned crystals 
of YBa2Cu307 and is roughly - 20 fiQ-cm. in the best untwinned crystals. A negative ex- 
trapolated value of Po means that the resistivity in the absence of superconductivity could 
not remain linear all the way down to zero temperature, but must turn over to some higher 
power law. The T^'s of twinned and untwinned crystals are about the same, despite the 
fact that the change in residual resistivity is of order 20 ^Q-cm. between the two types of 
samples. By contrast, if YBa2Cu307 were a d-wave superconductor, then, according to Fig. 
3, Tc would be strongly suppressed. 

This naive argument against the d-wave scenario is not strictly correct, since our calcula- 
tions have focused entirely on a single Cu02 plane, which implies that the residual resistivity 
in Fig. 3 represents only the planar contribution to po and does not include the contribution 
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from the CuO chains. For this reason, one can not simply take existing measurements of 
Tc and po in polycrystalhne or untwinned samples and infer unambiguous evidence for or 
against d-wave superconductivity; more systematic experiments are required. 

Currently, two methods could be used to perform such systematic experiments: substi- 
tution and irradiation. Substitutional studies [^,^,^ are able to introduce defects mainly 
in the Cu02 planes, as required, but find that the induced defects are generally magnetic. 
Abrikosov-Gor'kov-like behavior given by Eq. ([^) is then to be expected and provides no 
qualitative distinction between s- and d-wave pairing. Alternatively, defects can be induced 



by irradiation ||25| , |26|j27|l , but the location of these defects is often uncontrolled. Recent 
experiments p8|, though, suggest that low-energy electrons may preferentially disorder the 
Cu or O in the planes, which would enable a comparison of irradiation data with our pre- 
dictions. Current irradiation studies show that decreases linearly with fluence, as one 
would expect for weak disorder in the scaling curve of Eq. (|l|) , but it is not known how the 
fluence is correlated to the residual resistivity or whether the induced defects carry magnetic 
moments or not. Thus, although it is not yet possible to make a quantitative comparison 
of our results with current experiments, such a comparison should be possible with future 
systematic measurements. We hope that the present paper will serve as a stimulus for such 
experiments. 
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FIGURES 

FIG. 1. Critical temperature Tc normalized by the critical temperature in the absence of im- 
purities Tco vs. the bare impurity scattering rate T~^p in meV for the model of Ref. in the (a) 
strong-coupling and (b) weak-coupling cases calculated in the Fermi-surface-restricted Eliashberg 
formalism. The Tco's are 7.2 K (circle), 31.6 K (box), and 50.2 K (triangle); solid symbols denote 
strong-coupling results and empty symbols denote weak-coupling results. The solid lines are to 
guide the eye. 

FIG. 2. Normalized critical temperature T^/Tcq vs. pair-breaking parameter a = 
1/2(1 + Xz)TimpTco for the data in Fig. 1. Plot symbols are the same as in Fig. 1 with the 
solid symbols denoting strong-coupling calculations and the open symbols denoting weak- coupling 
calculations. For comparision, the analytic form of the scaling function [Eq. (jl])] is plotted as a 
plain solid line. Inset: Normalized critical temperature vs. pair-breaking parameter for the model 
of Refs. P] and [13| calculated from the exact Eliashberg equations with the coupling chosen so 



that Tco = 99.4 K. 

FIG. 3. Prediction of Tc as a function of the in-plane residual resistivity po due to non-magnetic 
impurities in a 90 K superconductor with a d-wave (shaded area) or an s-wave (dashed line) order 
parameter. The d-wave curve is computed from the generalized Abrikosov-Gor'kov form [Eqs. ([l|) 
and (^)] for experimental plasma frequencies oj*i between 1.1 and 1.4 eV, and the s-wave curve is 
simply a straight line due to Anderson's theorem. 
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